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Abstract 

Constitutive equations are developed for a polymer fluid, which is treated as a 
permanent network of strands bridged by junctions. The junctions are assumed to slide 
with respect to their reference positions under loading. Governing equations are derived 
by using the laws of thermodynamics under the assumption that the vorticity tensor for 
the flow of junctions is proportional to that for macro-deformation. Explicit expressions 
are developed for the steady elongational viscosity, as well as for the steady shear 
viscosity and normal stress functions. To verify the constitutive relations, three sets of 
experimental data are approximated on polystyrene solutions with various molecular 
weights. It is demonstrated that the model can correctly describe stress overshoot for 
the shear stress and first normal stress difference in start-up tests with various strain 
rates. Adjustable parameters in the governing equations change consistently with 
the strain rate, molecular weight and concentration of entanglements. To validate 
the constitutive equations, observations on low-density polyethylene melt in uniaxial 
extensional flow are compared with the results of numerical analysis when the material 
constants are found by matching experimental data in shear tests. 

Key-words: Polymer fluid, Non-affine network. Stress overshoot, Constitutive equations, 
Finite deformations 



1 Introduction 

This paper is concerned with modelling the rate-dependent response of polymer fluids at 
finite strains. In particular, the present study focuses on constitutive equations that can 
adequately describe stress overshoot under shear. The latter phenomenon has attracted 
substantial attention in the past three decades, because it provides a severe test for the 
analysis of rheological models. 
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Experimental studies demonstrate that under shear flow with a constant strain rate at 
isothermal conditions, the shear stress and the first normal stress difference grow mono- 
tonically with time to their ultimate values when the shear rate is relatively small. With 
an increase in the strain rate, the dependence of the shear stress on time becomes non- 
monotonic: it reaches a maximum at some instant (overshoot) and decreases to its limiting 
value afterwards (strain softening). At sufficiently high shear rates, the dependencies of both 
the shear stress and the first normal stress difference on time becomes non-monotonic, but 
the time when the first normal stress difference reaches its maximum noticeably exceeds that 
for the shear stress. The instants when the shear stress and the first normal stress difference 
become maximal are strongly affected by strain rate: the higher the shear rate is, the earlier 
the overshoots occur and the larger are the maximal stresses. 

Stress overshoot has been observed in conventional shear tests on polymer melts 121 El 
HlElini, immiscible molten blends of polymers [7j, polymer melts filled with micro- jH] and 
nano-particles I9i|10j, polymer solutions |lll |T2l CSI , solutions of associative copolymers [Hj, 
surfactant solutions [121, solutions of wormlike micelles jlSlEj; liquid crystalline polymers 
[TH] and discotic mesophases [THj, bicontinuous microemulsions [201, colloidal suspensions 
|2H I22j. and suspensions of rod-like macromolecules j2S|; to mention a few. This list demon- 
strates that stress overshoot is a rather common phenomenon observed in complex fluids. 
Comparison of experimental data for polymer solutions and melts with results of numer- 
ical simulation reveals that several constitutive equations can qualitatively describe stress 
overshoot under shear. We would mention among them the single integral model [2.., the 
Giesekus model j21I, the concept of temporary networks the finitely extensible nonlin- 
ear elastic (FENE) network model the reptation concept that incorporates segmental 
stretching [IIII2ZI, and the pom-pom model j2Hll2ni- 

However, even the most sophisticated constitutive equations that include a number of 
adjustable parameters can reproduce observations for the shear stress and the first normal 
stress difference in start-up shear tests only qualitatively, and, when their parameters are 
fitted to describe the stress overshoot, they fail to adequately predict the steady shear 
viscosity (see Figure 12 in [22]) ■ This observation may be explained by the fact that physical 
reasons for the origin of stress overshoot remain rather obscure, and not all of them are 
taken into account in conventional models. The following possible mechanisms for the onset 
of stress overshoot are mentioned: (i) changes in the concentration of entanglements and 
entanglement spacing jHj, (ii) changes in the orientation of polymer coils and nonuniform 
segmental stretching (iii) formation of an anisotropic mechanically-induced mesoscopic 
structure ^ that becomes unstable at relatively high strain rates JEj , and (iv) reorientation 
of meso-domains toward a steady-state aligned flow [IHl IHj . 

Although it is conventionally accepted that ordering of strands and their reorientation 
driven by shear flow provide the most important mechanism for stress overshoot and sub- 
sequent strain softening, to the best of our knowledge, no attempts have been made to 
incorporate this phenomenon into constitutive equations explicitly and to assess the effect of 
molecular structure of a polymer on the kinetics of reorientation. An important point here is 
that at finite strains, the evolution of an infinitesimal volume that can rotate with respect to 
macro- deformations is determined by two tensors: a symmetric rate-of-strain tensor and a 
skew-symmetric vorticity tensor (the so-called plastic spin), while the latter is conventionally 
disregarded in constitutive models for polymer fluids (despite the observation that "it is the 
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vorticity in shearing that leads to its extreme strain softening" 

The objective of this study is two-fold: (i) to develop constitutive equations for the 
response of polymer fluids that account for the evolution of the plastic spin, on the one 
hand, and that involve a relatively small number of material constants, on the other, and 
(ii) to demonstrate that these relations describe stress overshoot under shear of polymer 
solutions and melts quantitatively, while their adjustable parameters change consistently 
with strain rate and molecular structure of polymer fluids. 

A polymer fluid is treated as a network of chains bridged by junctions (entanglements 
and physical cross-links whose life-time exceeds the characteristic time of an experiment). 
The network is assumed to move non-affinely: the junctions between strands can slide with 
respect to their reference positions in the bulk material under deformation. Following [221 , we 
suppose that the vorticity tensor for viscoplastic flow of junctions is proportional to that for 
macro-deformation, whereas the rate-of-strain tensor for sliding of junctions is determined by 
the laws of thermodynamics. To simplify the analysis, we disregard rearrangement of strands 
in the network (which implies that we confine ourselves to relatively slow motions whose 
characteristic time exceeds the characteristic time for relaxation of stresses). This allows 
constitutive equations to be derived which describe any set of observations for the shear stress 
and the first normal stress difference in start-up shear tests by only four constants that have a 
transparent physical meaning. These parameters are found by matching experimental data 
for monodisperse polystyrene (PS) solutions with various mass-average molecular weights 
and molecular weights between entanglements and for a melt of low-density polyethylene 
(LDPE). 

The exposition is organized as follows. Basic kinematic relations for a non-affine network 
are derived in Section 2. Constitutive equations for a polymer fluid are developed in Section 3, 
where it is demonstrated that these relations generalize some well-known rheological models 
(Johnson-Segalman model [321; Phan-Thien-Tanner model jHSl; Giesekus model and 
Leonov model [SI])- Phenomenological relations for adjustable functions in the governing 
equations are suggested in Section 4. The governing relations are simplified for uniaxial 
extension in Section 5 and simple shear in Section 6, where explicit formulas are developed 
for the steady elongational and shear viscosities. In Section 7, material constants are found 
by fitting observations in start-up shear tests. A discussion of the effects of shear rate and 
molecular weight distribution on these quantities is provided in Section 8. The model is 
validated in Section 9 by comparison of results of numerical simulation for LDPE melt in 
extensional and shear fiows. Some concluding remarks are formulated in Section 10. 

2 Kinematic relations 

With reference to the concept of non-affine networks, a polymer fluid is modelled as a network 
of strands bridged by junctions that can slide with respect to their reference positions under 
deformation. Denote by F{t) the deformation gradient for transition from the initial state to 
the actual state at time t > 0, and by Fs(t) the deformation gradient for transition from the 
initial state to the reference state associated with sliding of junctions. It is presumed that 
the transformations described by the tensors F(t) and Fs{t) are volume-preserving, which 
means that the third principal invariants of these tensors equal unity. 
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The elastic deformation gradient Fe{t) [the deformation gradient for transition from the 
current reference state for junctions (the so-called intermediate state) to the actual state of 
the network] is determined by the conventional formula for the multiplicative decomposition 
of the deformation gradient 

F,{t)=F{t)-F:\t), (1) 

where the dot stands for inner product. Differentiation of Eq. (1) with respect to time t 
implies that 

r/F r/F 

-^{t) = L{t) ■ Fe(t) - F(t) ■ F:\t) • ^(t) • F;\t), (2) 



where 



(7F 

m = —{t)-F-\t) 



(3) 



is the velocity gradient. Using Eq. (1) and introducing the elastic velocity gradient Lc(t) 
and the velocity gradient for sliding of junctions ls{t) by the formulas analogous to Eq. (3), 



Le(t) 



dF, 



it)-F;\t), m 



dF. 



dt" " ^ ^' ^ dt 

we present Eq. (2) in the form 

Le(t) = L(t)-Fe(t)-l3(t)-F;l(t). 

Setting 

u{t) = F,{t)-ut)-F:\t) 

in Eq. (5), we arrive at the additive decomposition of the velocity gradients 

Le(i)=L(t)-L,(i). 
The left and right Cauchy-Green tensors for elastic deformation are given by 

Be(i) = Fe(i) • Flit), Ce(i) = Fj(i) • F,{t), 



(4) 

(5) 
(6) 

(7) 
(8) 



where T denotes transpose. We differentiate the first equality in Eqs. (8) with respect to 
time, use Eqs. (4) and (7), and find that 



dt 



it) 



L{t) ■ Be(i) + Be(i) • L^(i)J - [L,{t) ■ B,{t) + B,{t) ■ Lj{t) 

Differentiation of the other equality in Eqs. (8) results in 

dQ 



dt 



Kt)=2F:(i)-DeW-Fe(t), 



(9) 



(10) 



where 



De(t) = -[Le(t)+Lj(t) 



is the rate-of-strain tensor for elastic deformation. It follows from Eqs. (8) and (10) that 

dCz' 



dt 



-{t) = -2F;\t)--D,{t)-F;^{t). 



11) 



For volume-preserving deformations, the first and second principal invariants of the right 
Cauchy-Green tensor Ce{t) read 



Jei(t) = Xi(Ce(t)) = Ce(t) 1 1, Mt) = Ti{c;\t)) = C;\t) : I. 



(12) 



where Xi denotes the first invariant of a tensor, I is the unit tensor, and the colon stands for 
convolution. Differentiating the first equality in Eqs. (12) with respect to time and using 
Eqs. (8) and (10), we obtain 

^{t) = 2Be(t) : De(t). (13) 
Differentiation of the other equality in Eqs. (12) together with Eqs. (8) and (11) results in 

'-(t) = -2B;\t) : De(t). (14) 



dJe2, 



dt 

Equations (13) and (14) imply that the derivative of an arbitrary smooth function $( Jei, Je2) 
of the first two principal invariants of the right Cauchy-Green tensor for elastic deformation 
Ce{t) with respect to time t is given by 



'dt 



(Jei(t), Mt)) = 2[$i(t)Be(t) - Mt)^:\t)] : De(t), 



(15) 



where 



^m{t)^ j-—{jelit),Je2it)) (m=l,2). 



It follows from Eq. (7) that 
where 



De(t) = D(t)-D,(t), 



(16) 



D(t) = -[L(t) + L^(t) 



D,(t) = i[L,(t)+Lj(t) 



are the rate-of-strain tensors for macro-deformation and sliding of junctions, respectively. 
Combining Eqs. (15) and (16), we find that 



^{Mt),Mt)) = 2|[$i(t)Be(t)-$2(i)B,-^(t)] :D(t) 

'^i{t)B,{t)-Mt)B;\t)\ ■.B,{t)\. 



(17) 



It is convenient to present the velocity gradients L(t) and Ls(t) as the sums of the rate-of- 
strain tensors and vorticity tensors. 



L{t) = D(t) + n{t), L,{t) = u,{t) + n,{t), 



(18) 



where 



n,it)^-[U{t)-Lj{t) 



Without loss of generahty, we present the rate-of-strain tensor Ds(t) in the form 

B,{t) = (3{t)B{t) + M{t), (19) 
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where is a non- negative scalar function and 

M{t) = B,{t)-(3{t)B{t). (20) 
Substitution of expression (19) into Eq. (17) implies that 

= 2|(l - I3{t)) [$i(t)Be(t) - $2(t)B-i(t)]' : D(t) 

- $i(t)B,(t) - ^2mc\t)]' : M(t)}, (21) 



^(Jcl(t), Je2(t) 



where the prime stands for the deviatoric component of a tensor and we bear in mind that 
the tensors D(t) and M(t) are traceless. Following j32j, we assume the vorticity tensor for 
sliding of junctions to be proportional to the vorticity tensor for macro-deformation 

n,{t) =^it)n{t), (22) 

where 7(t) is a non-negative scalar function. Equations (18), (19) and (22) result in 

Ls(t) = M(t) + /3(t)D(t) + j{tMt). (23) 

Combining Eqs. (9) and (22), we arrive at the differential equation for the left Cauchy-Green 
tensor for elastic deformation 



it) 



L(t) ■ Be(t) + Be(t) ■ LT(t)J - Pit) [D(t) ■ Be(t) + Be(t) " D(t) 

M{t) ■ Be(t) + Be(t) ■ M(t)l - 7(t) \n{t) ■ Be(t) - B,{t) ■ n{t) 



where we take into account that M(t) is a symmetric tensor and f2(t) is a skew-symmetric 
tensor. Excluding the vorticity tensor Q{t) by means of Eq. (18), we obtain 



dB, 
~~dt 



it) = (l-7W)[L(t)-Be(t)+Be(t)-L^(t) 

+(7(t) - m) [d(^) ■ Be(t) + Be(t) ■ D(ty 

- M{t) ■ Be(t) + Be(t) ■ M(t) 



(24) 



Our aim now is to apply Eqs. (21) and (24) to derive constitutive equations for a polymer 
fluid by using the laws of thermodynamics. 



3 Constitutive equations 

Denote by the average number of strands per unit volume of a network and by w the 
average strain energy per strand. For an isotropic incompressible network, the quantity w is 
treated as a function of the first two principal invariants of the right Cauchy-Green tensor, 
w = w(^Jeiit), Je2(^))- Ncglcctiug the energy of interaction between strands (this energy is 
accounted for by means of the incompressibility condition j25]), we define the strain energy 
per unit volume as the sum of the strain energies of strands, 

W = Nw. (25) 
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For isothermal deformation of an incompressible medium at a reference temperature ©o; the 
Clausius-Duhem inequality reads 



Q{t) 



dW 



(t) + : D(t) > 0, 



where Q is internal dissipation per unit volume, and S is the Cauchy stress tensor. Substi- 
tuting expression (25) into this formula and using Eq. (21), we find that 



(26) 



Q{t) = -2N^(l-(3{t))[wi{t)B,{t)-W2{t)B;\t)\ : D(t) 

- Wi(t)B,(t) - W2{t)B;\t)]' : M(i)| + S'(f) : D(i) > 0, 

where the functions Wm{t) are given by 

dw 



(m = 1,2). 



(27) 



The dissipation inequality (26) is satisfied for an arbitrary deformation program, provided 
that the stress tensor reads 



S(t) = -P{t)l + 2iV(l - (3{t)) [wi{t)B,{t) - W2{t)B;\t) 

and the tensor M(i) is given by 

M{t) = a{t) [wi{t)B,{t) - W2{t)B;\t) 

Here P{t) is an unknown pressure, and a(t) is a non-negative scalar function. 
It follows from Eqs. (12) and (29) that 



(28) 



(29) 



1 



1 



M-Be + Be-M = 2a[wi(Bl- -J,^B,) -W2{l- -Je2B,) 

Substitution of this expression into Eq. (24) results in the kinetic equation 
dB 



dt 



■{t) = (l-7(t))[L(t)-Be(t) + Be(i)-L^(i)^ 

+ {^{t) - m) [D(^) • Be(i) + Be(i) • B{t) 

-2a{t)[w,{t)Bl{t) + ^{w2{t)J,2{t)-Wi{t)J,i{t))B,{t)-W2{t)l]. (30) 

Given a deformation program, Eqs. (28) and (30) provide a set of constitutive equations 
that involve four functions to be determined: the strain energy per strand w{Jci, Jc2), the 
function a(t) that characterizes the rate of internal dissipation, and the functions f3{t) and 
7(t) that establish connections between the rate-of-strain tensors and vorticity tensors for 
sliding of junctions and macro-deformation, respectively. It is worth noting that Eqs. (28) 
and (30) are apphcable for the description of the mechanical response of both solid polymers 
and polymer fluids, because no assumptions were introduced in their derivation regarding 
the specific properties of a network of strands. 
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An important feature of polymer fluids is that their time-dependent behavior reflects two 
different mechanisms at the micro-level: (i) sliding of junctions between chains with respect 
to their reference positions, and (ii) slippage of chains with respect to entanglements. Two 
approaches are conventionally employed to take into account the latter process (which is not 
included into the present model). 

According to the first (the concept of transient networks [23 ), if is assumed that strands 
are not permanently connected to their junctions, but detach from the junctions at random 
instants, see Figure 1 in When an end of an active strand separates from a junction, 
the strand is transformed into the dangling state, where the stress totally relaxes. When a 
free end of a dangling strand merges with a nearby junction at random time, the stress-free 
state of this strand coincides with the actual state of the network. 

According to the other approach (the reptation theory ^^), relaxation of stresses in 
chains occurs due to (i) their curvilinear diffusion along the tubes formed by surrounding 
macromolecules, and (ii) partial release from the tubes, see Figures 6.1 and 6.4 in 
As a result of this motion, the reference state of a strand differs from that of the vector 
that connects its junction points. This idea provides a physical ground for the pom-pom 
model where the average length of a strand in its stress-free state is introduced as an 
independent variable in the governing relations. 

Despite evident merits of these concepts, they share an important shortcoming: these 
approaches lead to (i) a substantial complication of the constitutive equations and (ii) a 
noticeable increase in the number of material functions and parameters. To avoid an un- 
desirable growth of experimental constants, on the one hand, and to make the constitutive 
equations thermodynamically consistent, on the other, we postulate that the characteristic 
rate of stress relaxation in strands (driven by their curvilinear diffusion in tubes or by their 
rearrangement in a transient network) substantially exceeds the rate of sliding of junctions. 
This implies that the number of strands per unit volume of a network should be treated 
not as a constant, but as a function of time whose values are rapidly (compared to the char- 
acteristic time of the deformation process) tuned to those that are determined by the current 
state of the network (an estimate for the relaxation time will be provided in Section 9). Our 
hypothesis is in accord with the double reptation concept jSH] that treat the response of a 
network of macromolecules as that induced by (i) rapid diffusive motion of chains in tubes 
(which is disregarded in the present model), and (ii) relatively slow relaxation of constrains 
imposed by surrounding macromolecules on the motion of a characteristic chain (which is 
accounted for in terms of the sliding process). 

The assumption regarding two relaxation processes with different time scales implies that 
the Clausius-Duhem inequality (26) remains valid for relatively slow deformation processes 
with A^ thought of as a function of time (according to the Doi-Edwards theory [33], N{t) 
is the number of strands per unit volume where stresses have not relaxed until instant t; 
in terms of the pom-pom model N[t) is proportional to the square of the stretch 
of a characteristic chain). According to this approach, the kinetic equation (30) remains 
unchanged, while the expression (28) for the stress tensor reads 



S(t) = -P{t)I + 2N{t)[l-P{t))[w,{t)B,{t)-W2{t)B^\t)\, (31) 
where N{t) is an adjustable function. 
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It is evident that the number of adjustable functions in Eqs. (30) and (31) is too large 
to find them with an acceptable level of accuracy by fitting experimental data in conven- 
tional rheological tests. To reduce this number, additional simplifications are conventionally 
introduced. 

1. Assuming strands to be Gaussian with the strain energy density 

w{Jcl, Jc2) = ^fJ'iJci - 3), (32) 

where /i is the average rigidity of a strand, setting a{t) = and •yit) = 0, and using Eq. 
(18), we arrive at the relations 



-P(t)I + iV(t)/iBe(t), 



(33) 



dt 



1 - f3{t)) [B{t) ■ Be(t) + Be(t) ■ D(t)J + [0(t) ■ Be(t) - Be(t) " n{t)\. (34) 

We now introduce the effective stress by the formula 

T(t) = iV(t)/iBe(t) - A(t)I, (35) 

where the function A{t) will be specified later. Differentiation of Eq. (35) with respect to 
time t results in 

dT dN , ^ ^ , ^ ^ dBc, ^ dk,^^ 

- = -W,B,W + N(t),^i,t) - -WI. 

Substitution of Eqs. (34) and (35) into this equality implies that 
dT' 

^(t) = {i-m)[m-m + m-m)\ + [m-m-m-m 



+^^WT(t) + 2(1 - m)mm) 



Kit) dN , , dA, ,1 , , 

The function A(t) is determined from the condition that the last term on the right-hand side 
of Eq. (36) vanishes, 

A(t) = CN{t), 
where C is an arbitrary constant. Setting 



1 dN 
N{t) dt 



it) 



1 



r(t)' 



where r(t) is a new material function (the characteristic time for rearrangement of tubes 
according to the reptation theory), we find that 

N{t) = Ar(0)exp[- f\{s)ds 



Substitution of these expressions into Eq. (36) results in the constitutive equation of the 
Johnson-Segalman model |32j . 



dT 

— it) = (l-/?(t))[D(t)-T(t) + T(t)-D(t) 
1 



+ 



n{t) ■ T{t) - T(t) ■ n{t) 



'r{t) 



2r/(t)D(t) - T(t) 



(37) 



where the viscosity rj is given by 

rj{t) = CNq (l - r (t) exp [- ^* T{s)ds 
2. Adopting hypothesis (32) and setting a{t) = and 7(t) =0, we find that 

E(t) = -P(t)I + iV(t)/i(l-/3(t))Be(t), 
Bf(t) = -/5(t)[D(t)-Be(t) + Be(t)-D(t)], 

where 



Bf(t) 



rfB, 
~dt 



(t)-L(t)-Be(t)-Be(t)-L^(t) 



is the Oldroyd contravariant derivative. Introducing the effective stress T(t) by Eq. (35), 
we obtain 



T^(t) 



-/?(t)[D(t)-T(t)+T(t)-D(t) 
-2(l-/3(t))A(t)D(t) + 



1 dN , , 

+ T77TT — (t)T(t) 



iV(t) rft 



I. 



Equating the last term in the right-hand side of this formula to zero, we arrive at the Phan- 
Thien-Tanner constitutive model 



T^(t) = -/5(t) D(t) ■ T(t) + T(t) ■ D(t) 



1 



Tit) 



2r^(t)D(t) - T(t) 



(38) 



3. Accepting formula Eq. (32) for the strain energy of a Gaussian strand and assuming 
that [3{t) = and ^{t) = 0, we find that the Cauchy stress tensor S(t) is given by Eq. (33), 
where the function Be(t) is governed by the differential equation 



B^it) = -a{t)fx B'iit) - ^Jei(t)Be(t) 



Substitution of expression (35) for the effective stress T{t) into this equahty implies that 



+2A(t)D(t) + 



rA(t) dN . . a(t) 1 

; , (t) - ttttA^ t - -a(t)uA(t)J^i(t) 
N{t) dt ^ ' N{t) ^ ' 3 V ;p V ; cu ; 



dk 



it) 



I. 



Equating the last term of the right-hand side of this equation to zero, which results in the 
differential equation for the function A{t), 



dA , ^ . / N 



1 dN 



it) 



-N{t) dt 

and introducing the new adjustable functions 



^A(t) - ia(t)/iJei(t) 



-it) 



1 dN , , a(t)A(t) 1 ^ , , 



7]{t)=A{t)r{t), 
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we arrive at the formula 



^ ' N{t) ^ ' 



r{t) 



2r^(t)D(t) - T(t) 



(39) 



which coincides with the Giesekus constitutive equation 

It is worth noting that the material functions r(t) and ri{t) in Eqs. (37) to (39) are not 
arbitrary, but they obey some additional restrictions. This is caused by the fact that we 
consider relaxation of stresses in strands driven by their diffusion inside tubes as a rapid 
process compared to rearrangement of tubes (treated as sliding of junctions). An advantage 
of our approach is that it operates with the Cauchy deformation tensors conventionally 
employed in continuum mechanics (no configurational tensors are introduced) and is based 
on the classical multiplicative decomposition of the deformation gradient (unlike constitutive 
equations that explicitly use a slip tensor). 

4. Finally, under the assumptions that (i) w is a symmetric function of its arguments, 
which implies that 

Wi{t) = W2{t) = W{t), 

and (ii) the functions P(t) and •yit) vanish, Eqs. (28) and (31) are transformed into the 
Leonov constitutive model 



m 
Kit) 



-Pit)I + 2N{t)w{t)(B,{t)-B-\t)] 
-2a{t)w{t)[Bl{t) + i(je2(t) - Jcl(t))Be(t) 



I 



Unlike the above approaches, this study focuses on the case when the functions l3{t) and 
7(t) coincide, 

y{t)=P{t)=bit). (40) 

This allows us to omit the second term on the right-hand side of Eq. (30), which may result 
in stress oscillations at simple shear that have no physical meaning. From the kinematic 
standpoint, condition (40) means that the velocity gradients for sliding of junctions and 
macro- deformation are connected by the linear equation, see Eqs. (18) and (23), 

U{t) = b{t)L{t)+M{t), 

where M(t) is a symmetric traceless tensor given by Eq. (29). Substitution of expression 
(40) into Eq. (30) implies the differential equation 



~dt 



it) 



l-6(t))[L(t)-Be(t) + Be(t)-LT(t) 

-2a{t)[wi{t)Bl{t) + ^[w2it)J,2it) - Wi{t)J,^{t))B,{t) - W2it)l 



(41) 



For a network of Gaussian strands, see Eq. (32), the stress-strain relations (31) and (41) 
read 



dB^ 
~dt 



it) 



-P{t)I + G{t)(l-b{t))B,{t), 
l-6(t))[L(t)-Be(t) + Be(t)-LT(t) 



a{t) Bl{t) - -Jei(t)Be(t) (42) 
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with 

a{t) = a{t)fx, G{t) = fxN{t). 

Our aim now is to specify the dependencies of a, b and G in Eqs. (42) on parameters that 
characterize macro-deformation of a polymer fluid. 



4 Adjustable functions 

1 

We introduce the rate intensity Di by the conventional formula Di — ^2D : ^ and adopt 
the following equation for the elastic modulus G: 

^""^ exp{KGD^-), (43) 



In Go 

where Gq, Kq and are non-negative material constants. According to Eq. (43), Gq is the 
plateau modulus at small strain rates, whereas Kq and vq characterize a strong increase in 
G with strain-rate intensity. 

The parameter h is assumed to obey the power law, 

b = K^D^" (44) 

with non-negative adjustable parameters and Uh. Formula (44) means that no sliding of 
junctions occurs at very small strain rates, and the rate of sliding grows with the rate of 
macro-strain. Equations (43) and (44) are applicable within a limited range of strain rates, 
as the modulus G and the parameter b should reach some limiting values at sufficiently high 
strain rates. 

The coefficient a is split into the product of two functions 

a = C^, (45) 

where the function ^ characterizes the effect of strain and the function A determines the 
effect of strain rate on the rate of energy dissipation [described by means of the coefficient 
a that is proportional to a in Eqs. (26) and (29)]. 

The dependence of the parameter A on strain-rate intensity is determined by the 
power-law similar to Eq. (44), 

A = Ao + KaD'^^, (46) 

where Aq, Ka and ai'e non-negative coefficients. Eq. (46) means that A monotonically 
increases with strain rate from its plateau value Aq. Unlike Eq. (44) for the function b{D[), 
we do not presume A to vanish at small strain rates; a reason for this claim will be explained 
later in Section 9. 

It seems natural to suppose that for a developed flow with large strains, the motion of 
junctions at the micro-level is affine with macro-deformation. It follows from Eqs. (23) and 
(40) that this hypothesis is tantamount to the assertion that the coefficient a in Eq. (29) 
vanishes. This conclusion implies that C should approach zero at large deformations (that 
is at large strain energies). To approximate such a dependence, the exponential function is 
chosen, 

C = exp[-r(Jei-3)l, (47) 
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where r > is a constant. Our choice of the term J^i — 3 in Eq. (47) may be justified by 
the fact that this expression is proportional to the strain energy of a Gaussian strand (32). 
Combining Eqs. (45) to (47), we arrive at the formula 

a^{Ao + KaD^^) exp [-r ( Jei - 3)] . (48) 

Formally, Eqs. (42) to (44) and (48) involve 9 adjustable parameters: Aq, Gq, Ka, Kf,, Kq, 
ua, I'h-i and r. It should be noted, however, that any set of experimental data at an 
arbitrary deformation with a constant strain-rate intensity Di is determined by 4 constants 
only: A, b, G and r. The other quantities can be found by fitting the functions A{Di), 6(A) 
and G{Di) by Eqs. (43), (44) and (46). 

Our aim now is to apply constitutive equations (42) for the analysis of uniaxial extension 
and simple shear of a polymer fiuid. 



5 Uniaxial extension 

Uniaxial extension of an incompressible medium is described by the formulas 

xi = k{t)Xi, X2 = k~^t)X2, xs = k~^{t)X3, 



(49) 



where {Xi} are Cartesian coordinates in the initial state, {xi} are Cartesian coordinates in 
the deformed state, and k = k{t) stands for elongation. We suppose that transition from the 
initial state to the current reference state that determines the sliding process is described by 
the equations similar to Eqs. (49), 



6 = «:(t)Xi, ^2 = f^'Ht)X2, ^, = n-Ht)X3, 



(50) 



where {^i} are Cartesian coordinates in the intermediate state, and k, = K,{t) is a function to 
be found. According to Eqs. (49) and (50), the deformation gradients F and Fg are given 

by 

F^keiBi + k 2(6262 + 6363), Fs = K6i6i + K 2(6262 + 6363), (51) 
where are unit vectors of the frame {Xi}. It follows from Eqs. (1), (8) and (51) that 



/k\2 K. , 

(- 6161 + -(6262 + 6363) 



and 



J. 



cl 



k\2 

k' k 



The velocity gradient L, the rate-of-strain tensor D, and the strain-rate intensity D\ read 



D = - 1^6161 - -(6262 + 6363^ 



\k\ 
k 



73, 



(52) 



where k — dk/dt. STibstitution of these expressions into Eqs. (42) results in the differential 
equation for the function K{t), 



K 
K 



, k a 



(53) 
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The Cauchy stress tensor S is given by 

S = Sieiei + 112(6262 + 6363), 
where the stress difference AE = Ei — E2 reads 

AE = G(l-6)[(^)^-^]. (54) 

Given a deformation program k{t), Eqs. (53) and (54) together with the phenomenological 
equations (43), (44) and (48) determine the stress difference AE as a function of time. 
For uniaxial extension with a constant rate of Hencky strain 

k 

1 = '^ 

it is natural to search a solution of Eq. (53) in the form 

K(i) = ^(t)exp(et), (55) 

where z(t) is a function to be found. Substituting expression (55) into Eq. (53) and using 
Eq. (48), we find that 

i + (1 _ = ^(1 - z') exp(-ri^:^4^) ' ^(0) = 1' (^6) 

6z ^ z^ ' 

where A, h and G depends on e. It follows from Eq. (56) that when t — > cxd, the function 
z{i) approaches its limiting value zq^ which is determined by the nonlinear equation 

iz^ oxp(-ri^^|±M) = (57) 
According to Eqs. (54) and (55), the transient elongational viscosity 

rf. - ^ (58) 

is given by 

tz^ 

Equation Eq. (57) implies that when i ^ 00, the transient viscosity r]^ tends to its limiting 
value 

3G, 2 ( 1 — 3^0 + 22^0 \ 
77e = ^(1 - hf exp(r ^ 

This formula together with Eq. (57) determine the steady elongational viscosity 
function of the Hencky strain rate e. When r = (no effect of strain on the coefficient a), 
we find that 

on 

^^ = —{\-h)\ (59) 

According to Eqs. (44) and (48), Eq. (59) describes a pronounced decrease in the steady 
elongational viscosity at large strain rates observed in numerous experiments. 
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6 Simple shear 

Simple shear of an incompressible medium is described by the equations 

Xi ^ Xi + k{t)X2, X2^X2, X3^X3, (60) 

where {X^} and {xi} are Cartesian coordinates in the initial and actual states, respec- 
tively, and k{t) stands for shear. We describe transition from the initial to the intermediate 
configuration that characterizes sliding of junctions as a superposition of simple shear and 
three-dimensional extension, 

Cl^Xl{t)Xi + K{t)X2, 6 = A2(i)X2, C3^Ht)X3, (61) 

where {^j} are Cartesian coordinates in the stress-free state, and \i{t) and K{t) are functions 
to be found. The functions Xi{t) obey the incompressibility condition 

A1A2A3 = 1. 

It follows from Eqs. (60) and (61) that the deformation gradients F and Fg read 

F = eiei + 6262 + 6363 + A;eie2, 
Fg = Aieiei + A2e2e2 + Xsese^ + neie2. 

Substitution of these expressions into Eq. (1) implies that 

Fe = pieiei + ^26262 + ^36363 + 06162, (62) 

where 

Pi = A2A3, p2 = AiA3, P3 = AiA2, (f)^ XsiXik - k). 
In the new notation, the incompressibility condition is given by 

P1P2P3 = 1- (63) 

It follows from Eqs. (8) and (62) that the left Cauchy-Green tensor for elastic deformation 
is determined as 

Be = (Pi + 0^)6161 + P26262 + P36363 + P20(6l62 + 6261) . (64) 

The velocity gradient L, the rate-of-strain tensor D, and the strain-rate intensity A read 

1 • 

L = A;6i62, D = -A;(eie2 + 6261), A = k. (65) 

We substitute expressions (64) and (65) into the second equation in Eqs. (42). Omitting 
simple but tedious algebra, we arrive at the following equations: 

Pi = -^Pi(2p? -P2 -Ps - 0^)> Pi(0) = 1, 
6 

P2 = ~P2{-pi + 2pl-pl + 2<f>^), P2(0) = l, 

h = ~Pli-pl-pl + 2pl-(p^), P3(0) = l, 

= {l-b)kp2-^<P{5pl + 2pl-pl + 24>'), 0(0) = 0. (66) 
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It is easy to check that Eqs. (66) imply the incompressibihty condition (63) for an arbitrary 
deformation program k{t). Combining Eqs. (42) and (64), we obtain 



S = SiieiCi + 2226262 + 2336363 + Si20(6ie2 + 6261), 



where 



S12 



= -P + G{l-b)pl, 
= -P + G{l-b)pl 
= (7(1-6)^20. 



It follows from these equations that the first normal stress difference - 
the second normal stress difference A^2'^ = S22 — S33 read 

7V+ = G{1 - bM -pI + <f>'), 7V+ = G(i - b){pl - pI). 



(67) 
S22 and 

(68) 



Equations (66) to (68) describe the time-dependent response of a polymer fluid under shear 
with an arbitrary deformation program k{t). Our aim now is to analyze these equations for 
shearing with a constant strain rate k and to study their steady-state solutions at t ^ 00. 
Assuming that Pi{t) piQ and (f){t) — > 0o as t — > 00, we find from Eqs. (66) that the limiting 
values pio and 0o satisfy the algebraic equations 



2 _ 2 

P20 PSQ 



= 0, 

-P?o + 2P2O-P3o + 20O = 0, 

„2 2 I o„2 i2 



-Pio - P20 + 2P30 - (Po 



0. 



(69) 



To derive these equations, the incompressibihty condition (63) is employed which excludes 
the case when any of pio vanishes. Subtracting the third equality in Eqs. (69) from the first, 
we find that 

(70) 

(71) 



2 2 

PSQ — PlO- 



Substitution of expression (70) into the first equahty in Eqs. (69) results in 



pIo 



pIo + 'Po- 



Combining Eqs. (63), (70) and (71), we obtain 

P2o{pIo + (Pl) = 1- 

It follows from Eqs. (66) that 4>q satisfies the equation 

X 

(1 - b)kp2o = -Co^0o(5pio + 2p2o - pIo + 
where, according to Eqs. (46) and (64), 



(72) 



Co 



exp 



-r{plo+plo+Plo + <Pl-^) 
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Substitution of expressions (70) and (71) into these equations implies that 

{l-b)k 



MP2o + (l^o) 



-P20, 



Co = exp[-3r(^p2o + 



- 1 



Excluding the sum + 4'o by means of Eq. (72), we find that 



{l-b )k 
1 



-P20 exp 



3r{p. 



20 
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(73) 



Finally, substitution of expression (73) into Eq. (72) results in the transcendental equation 
for P20, 



1 + 



-P2oexp(6r(p2o^ - 1] 



1. 



(74) 



It follows from Eqs. (67) that the shear stress S12 tends to its limiting value 

S120 = ^(1 - b)p2o4>o 

when t approaches infinity. Substituting expression (73) into this formula and introducing 
the steady shear viscosity 

S12 



we obtain 



k 



G ■ V 

Vs = ^(1 - bfploexp[3r{p-^ - 1) 



(75) 



In the case when r = 0, which is equivalent to the neglect of the infiuence of strain on the 
coefficient a, Eq. (75) reads 

Vs = j{l-brpl,. (76) 
Comparing Eqs. (59) and (76) and introducing the Trouton ratio 

Tr = — , 

we obtain 

Tr = 3p2o'- 

Combining this formula with Eq. (74) (where we set r = 0) and bearing in mind Eq. (65), 
we arrive at the equation for the Trouton ratio 



Tr = 3 1 



r3i(l -6)A 



A 



Tr~3 



This formula implies that under the condition < 1, the Trouton ratio equals 3 at very 
small strain rates, and it increases with the strain-rate intensity, at least at relatively small 
Di, in agreement with the experimental data reported in |37j . 

We return now to the general case r > and calculate the limiting values Ni and N2 of 
the normal stress differences when t 00. It follows from Eqs. (68), (70) and (71) that 



A^i = 2G(1 - 6)0; 



0' 



iV2 = -G(l-6)02. 
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Introducing the steady normal stress functions = N^k ^ and using Eq. (73), we find that 

2C 1 

^1 = ^(1 - hfplo^M^rip^o - 1)1, ^2 = --^1. (77) 
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Equations (74), (75) and (77) together with the phenomenological relations (43), (44) and 
(46) determine the steady shear viscosity and normal stress functions for an arbitrary shear 
rate k. Equation (77) correctly predicts that the second normal stress coefficient \E'2 is 
negative and proportional to the first normal stress function It is conventionally assumed 
for polymer fiuids that the coefficient of proportionality between and '^2 should be of order 
of 0.1 to 0.3 |37i, but the ratio |\E'2|/\E'i is strongly affected by experimental conditions, and 
its precise value remains unknown jHE]- Recent experiments on suspension in viscous fluids 
show that this ratio is located in the interval between 0.1 and 0.7 with the most probable 
value of 0.5 for highly viscous liquids [221 • The latter result is in excellent agreement with 
Eqs. (77). 



7 Fitting of observations 

Our aim now is to demonstrate that constitutive equations (42) can correctly approximate 
observations in start-up shear tests with various strain rates and to show that the adjustable 
parameters change consistently with shear rate and characteristics of molecular weight dis- 
tribution. For this purpose, we focus on three sets of experimental data for solutions of 
polystyrene in tricresyl phosphate. For a detailed description of the experimental procedure 
and the material properties, we refer to the original studies IT^ IT!^ . The main parame- 
ters of the solutions (abbreviated here as PSl, PS2 and PS3) are collected in Table 1. The 
polydispersity index of PS2 and PS3 are not provided, but it is mentioned that they have 
"sharp molecular weight distributions" OIISI- Based on these observations, we treat the 
solutions as monodisperse and concentrate on the effects of mass-average molecular weight 
and the number of entanglements per chain M^jM^ [M^ stands for the average molecular 
weight between entanglements) on material constants. 

We begin with fltting the observations on PSl for the shear stress S12 and the flrst normal 
stress difference depicted in Figures 1 and 2. These figures show that the functions 
Ei2(t) and Ni~{t) increase monotonically with time at small shear rates k and demonstrate 
pronounced stress overshoots at relatively large k. As the experimental data were reported 
in Figure 4 of [llj with the use of the regular time-scale, we follows the same approach. 
The other sets of data (see Figures 6 to 8 below) are presented in the logarithmic time-scale 
(log = logj^p) that makes more transparent the characteristic features of stress overshoot. 

First, we approximate the experimental dependencies Si2(t) and Ni~{t) measured at the 
highest strain rate fcmax = 4.0 s~^, by using four constants. A, b, G and r [the coefficient a in 
Eqs. (66) is calculated from Eq. (45), where ( is given by Eq. (47)]. To find the coefficients 
A, b, G and r, we fix some intervals [0,y4inax], [0,&max], [0,G'max] and [0,rmax], where the 
"best-fit" parameters A, b, G and r are assumed to be located, and divide these intervals 
into J subintervals by the points A^'^^ = iAA, b^^^ = jAb, G''''^ = kAG and r^^^ = lAr 
{i,j,k,l = 1,. . . ,J- 1) with AA = A^ax/ J, Ab = b^^J J, AG = G^^J J and Ar = r^ax/ J- 
For any set {y4^*\ fe^-'^ r*^^)}, Eqs. (66) are integrated numerically by the Runge-Kutta 
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method with the time-step At = 1.0 ■ 10"'^ and the quantities S12 and are calculated 
from Eqs. (67) and (68). The "best-fit" parameters A, b, G and r are determined from the 
condition of minimum of the function 



R 



E 



on the set {A^'^\ h^^\ G^^\ r*^') {i,j,k,l = 1, . . . , J — 1)}. Here the coefficient 5 characterizes 
a "weight" of the second term compared to the first one, the sum is calculated over all 
experimental points tm depicted in Figures 1 and 2, Ei2^ and Ni^"^ are the shear stress 
and the first normal stress difference measured in a test, and E™™ and N^^^^ are given 
by Eqs. (67) and (68). In the numerical analysis, we set J = 10 and S = 0.05. After 
finding the "best-fit" values A^^\ h^^\ G^''^ and r*^'\ this procedure is repeated twice for the 
new intervals [A^'-^\ A^'+% [b^^-^\b^^+% [G^'^-i), G^'^+i)] and [r'-i), r('+i)], to ensure an 
acceptable accuracy of fitting. 

After finding the "best-fit" value of r by matching the data obtained at fcmax, we fix this 
parameter and approximate the observations at other shear rates with the help of only three 
material constants. A, b and G. Each pair of curves (for S12 and Ni) measured at a given 
shear rate k is matched separately. 

Figures 1 and 2 demonstrate good agreement between the experimental data and the 
results of numerical simulation. The model correctly predicts that the shear stress E12 
reaches its maximum before the first normal stress difference approaches its highest 
value. It is worth mentioning that the quality of fitting observations in Figures 1 and 2 is 
substantially better than that reported for other constitutive models in [221 (underprediction 
of the shear stress) and in [27j (overprediction of S12). 

The quantities A, b and G are plotted in Figures 3 to 5 versus shear rate k. The function 
A{k) in Figure 3 is approximated by the simplified version of Eq. (46) with Aq = 0. Bearing 
in mind expression (65) for D^, we present Eqs. (44) and (46) in the form 



log A = log Ka + T^A log k, log b = log + Ub log k. 



(79) 



The coefficients Ka, Kb, va and Vb in Eqs. (79) are determined by the least-squares technique. 
The function G{k) in Figure 5 is fitted by Eq. (43), where the coefficients Gq, Kq and are 
found by the nonlinear regression algorithm. These quantities are listed in Table 2. Figures 
3 to 5 show reasonable agreement between the experimental data and their approximations 
by Eqs. (43) and (79). 

We proceed with fitting observations for polystyrene solutions PS2 and PS3, where we 
confine ourselves to matching experimental data for the shear stress S12 only. The same 
algorithm of fitting is employed as for PSl, with the only exception that the last term in 
Eq. (78) is omitted (5 = 0). 

The experimental dependencies Si2(t) measured at various strain rates k are plotted 
together with the results of numerical simulation in Figure 6 for PS2 and in Figure 7 for 
PS3. These figures demonstrate good agreement between the observations and the results 
of numerical analysis at all strain rates under consideration. It is worth noting that the 
model can correctly reproduce not only stress overshoot, but also weak stress undershoot 
revealed in experiments. The adjustable parameters A, b and G are presented in Figures 3 
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to 5 as functions of shear rate together with their approximations by Eqs. (43) and (79). 
The "best- fit" values of the parameter r are collected in Table 2. 

Our aim now is to discuss the effects of strain rate and molecular weight of polystyrene 
solutions on the adjustable parameters in the constitutive equations. 

8 Discussion 

According to Figure 5, the elastic modulus G grows with strain-rate intensity A for all 
solutions. It follows from Table 2 that the rate of increase in G [the exponent in Eq- (43)] 
monotonically decreases with number of entanglements per chain M^jMf,. This observation 
appears to be natural if we treat changes in G with strain rate within the Doi-Edwards 
reptation theory as a result of rearrangement of tubes driven by macro-deformation. The 
latter implies that changes in the number of non-relaxed strands with k should be extremely 
pronounced for a weakly entangled network (where a strand is practically the same as a 
chain), and should be relatively weak for a strongly entangled system (where the number of 
entanglements per chain is rather large, and only strands located at the ends of a chain can 
easily escape from the tubes). 

The other parameter that characterizes the effect of strain rate on the elastic modulus, 
the coefficient Kq in Eq. (43), appears to correlate with mass-average molecular weight M^: 
the higher the molecular weight is, the smaller is the pre-factor Kq. This implies that both 
factors, the mass-average molecular weight and the average number of entanglements per 
chain strongly affect the dependence of the shear modulus on strain rate. 

Table 2 shows no direct correlations between the plateau modulus Go and the mass- 
average molecular weight M„. It reveals, however, a noticeable dependence of Gq on the 
ratio M^/Mc- the larger the number of strands per chain is, the higher is the modulus Gq. 
This conclusion is in excellent agreement with the classical theory of rubber elasticity, which 
predicts that the elastic modulus is determined by the number of strands per unit volume, 
and it is independent of the average length of a strand. 

According to Figure 3, the parameter A grows with strain-rate intensity, and the rate 
of increase in A with shear rate is practically independent of the network structure (the 
exponent va ~ 0.4 weakly depends on molecular weight and concentration of entanglements). 
Table 2 demonstrates, however, that the absolute values of A are noticeably affected by the 
ratio M^/Me- the coefficient Ka substantially decreases (by an order of magnitude) with 
average number of entanglements per chain. This finding seems natural if we recall that A 
(which is proportional to a) characterizes the influence of the network elasticity on the non- 
affine flow of junctions that is described by the tensor M. When the number of junctions 
per unit volume is relatively large, stresses in strands are not sufficient to substantially 
affect the sliding process (which implies that Dg ~ /3D). On the contrary, for a network 
with a relatively small concentration of junctions, stresses in strands strongly influence flow 
of junctions (which means that the tensor M does not vanish, and, as a consequence, the 
coefficient a in Eq. (29) is relatively large). 

According to Figure 4, the parameter 6, that characterizes the rate of sliding of junctions 
driven by macro-deformation, grows with shear rate at all strain rates under consideration. 
For all solutions, the parameter h remains below unity. The rate of increase in h with k 
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reaches maximum for PS2 (the lowest M„), whereas these rates are rather modest for PSl 
and PS3 (solutions with high molecular weights). Table 2 shows that the exponent Uh in Eq. 
(44) increases with mass- aver age molecular weight, and it reveals no correlations with the 
average number of entanglements per chain. 

According to Table 2, the coefficient r noticeably decreases with the ratio M^/M^,. This 
implies that the effect of strain energy on the rate of internal dissipation [characterized by 
the coefficient a in Eq. (29)] is substantial for solutions of weakly entangled polymers and 
may be disregarded for strongly entangled networks. 

9 Validation of the model 

To verify the constitutive equations, we apply two approaches. According to the first, we 
analyze the same deformation mode (shear flow), but for a noticeably denser system and 
compare the conclusions drawn in the previous section with results of numerical simula- 
tion. Following the other approach, we focus on the response of a polymer melt at another 
deformation mode (extensional flow) and compare predictions of the governing equations 
(with the adjustable parameters found by fltting experimental data under simple shear) 
with observations. 

1. We begin with the flrst approach and formulate explicitly some assertions that can 
be drawn for a dense polymer system from the analysis of results reported in Figures 3 to 5 
and Table 2. As the network under consideration we chose LDPE melt (lUPAC A) whose 
physical properties are summarized in Table 1 with reference to jlOl (the mass-average 
molecular weight and the number- average molecular weight Mn of the melt are given in 
jlU] : to assess the average number of junctions along a chain, we use the graph depicted in 
Figure 3 of jH] and presume that the number of junctions exceeds the number of branching 
points for the highly branched LDPE melt). Our choice may be explained by two reasons. 
First, the observations in start-up shear tests on LDPE melt at the temperature Qq = 150 °C 
reported in are accompanying by those in transient extensional tests on the same polymer 
at the same temperature [32]; which allows results of numerical simulation to be compared 
at two different deformation modes. Secondly, these observations have recently been used to 
validate the pom-pom model |2H] • Thus, the quality of fltting by the constitutive equations 
(42) and the pom-pom model may be collated on the same set of data. 

The following hypotheses are introduced regarding the adjustable parameters of the melt: 

(i) The modulus G should substantially exceed that for the solutions at all strain rates 
(this fact is not trivial because G is determined by fltting the entire dependencies of 
the shear stress on time). The coefficient Go should be relatively large, whereas the 
exponent ug in Eq. (43) should be rather small (these quantities are determined by 
the ratio M^/M^ that is high for the melt). The pre-factor Kg in Eq. (43) should be 
relatively large (this quantity is inversely proportional to the mass-average molecular 
weight, while M„ of the melt is small compared with the molecular weights of PS 
solutions). 

(ii) The coefficient r in Eq. (47) that describes the effect of strain energy on the rate of 
sliding of junctions should be relatively small (because the ratio M^/M^ of the melt 



21 



exceeds that of PS solutions). 

(iii) The exponent ua of the meh should be close to those of PS solutions (this parameter is 
independent of the molecular structure). The coefficient Ka in Eq. (46) that accounts 
for the influence of strain rate on the rate of sliding of junctions should be quite small 
(this quantity is inversely proportional to the average number of strands in a chain). 

(iv) The exponent z/j, in Eq. (44) that describes the effect of strain rate on the parameter b 
(the latter characterizes the proportionality between the velocity gradients for macro- 
deformation and sliding of junctions) should be relatively small (because the ratio 
Mw/Me of the melt exceeds that of PS solutions). 

To examine these hypotheses, we begin with the approximation of the experimental data 
for the shear stress S12 plotted as a function of time in Figure 8. To reduce the number of 
material constants, we set r = in accord with assumption (ii). To find the coefficients A, 
b and G, each curve Si2(t) is fitted separately by using the same approach that was utilized 
in Section 7. We chose some values of A, b and G from the intervals [0,ylinax], [0, &max] 
and [0, Gmax], where the "best-fit" parameters are assumed to be located, integrate Eqs. 
(66) numerically, and calculate the shear stress E12 by means of Eq. (67). The "best-fit" 
coefficients A, b and G are found from the condition of minimum of the cost function R in 
Eq. (78), where the last term is disregarded. 

Figure 8 demonstrates excellent agreement between the observations and the results of 
numerical simulation. It is worth noting that the quality of fitting the experimental data 
is noticeably higher than that for the pom-pom model with a substantially larger number 
of experimental constants, compare Figure 8 with Figure 9 in |28j. By no means, this 
conclusion imphes that Eqs. (42) are superior with regard to the constitutive equations in 
the pom-pom model, as the latter are grounded on the solid basis of the molecular theory 
of polymer networks. We suppose that the good quality of matching observations with the 
smaller number of adjustable parameters revealed in Figure 8 may be explained by the fact 
that a detailed description of orientation and stretch of chains is excessive for the analysis 
of shear fiows with relatively large strain rates. 

The modulus G determined in the fitting procedure is plotted versus strain rate k in 
Figure 5. The dependence G{k) is matched by Eq. (43) with the adjustable parameters 
Go, Kg and reported in Table 2. Comparison of the results hsted in Table 2 shows 
that the modulus G of the melt substantially exceeds that for the solutions (at least, by 
three orders of magnitude), the exponent for the melt is small compared to that for PS 
solutions, whereas the pre-factor Kq exceeds that for the solutions. All these conclusions 
are in perfect agreement with assumption (i). 

It is obvious that assumption (ii) is fulfilled, as the value r = was set in the fitting 
procedure. 

The coefficient A is plotted versus k in Figure 3 together with its approximation by Eq. 
(79) with the coefficients Ka and ua listed in Table 2. Comparison of these parameters with 
appropriate parameters for PS solutions demonstrates that assumption (iii) is satisfied. 

The parameter b is plotted versus strain-rate intensity D[ in Figure 4. This figure demon- 
strates that b is practically constant, which means that the exponent z/f, in Eq. (44) vanishes, 
in agreement with assumption (iv). 
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This analysis reveals that the approximation of experimental data in shear tests for 
polymer solutions by the constitutive model implies the assumptions regarding the behavior 
of material parameters of a melt that appear to be physically plausible. 

2. To further examine the constitutive equations, we intend to deduce a relationship that 
does not contain adjustable parameters and that connects some quantities which may be 
measured directly in conventional tests. For this purpose, we confine ourselves to relatively 
dense polymer systems (melts), and, based on assertion (ii), set r = in the governing 
equations. According to Eqs. (65) and (77), in this case the first normal stress function 
reads 

2G{l-hf 



^1 



0) 



Equations (76) and (80) result in 
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^1 _ 2 

On the other hand, Eqs. (65) and (74) (where we set r = 0) imply that 



Combining these equalities, we find that 

^ 2r^ J ' 
2Vs 



P20 



which means that 
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It follows from Eqs. (59) and (76) that 



rje = 3r/sP2o ■ 

Substituting expression (81) into this equality, we arrive at the formula 



Ve = 3?7s 



^ 2r]s ' 



12) 



which expresses the steady elongational viscosity rje in terms of the steady shear viscosity ?7s 
and first normal stress function \E'i for an arbitrary strain-rate intensity D^. 

To check the validity of Eq. (82), we use observations on polypropylene melt (M„ = 
3.5 ■ 10^ g/mol) at Gq = 190 °C reported in Figure 10 of It should be noted that 

steady extensional and shear tests are conventionally performed at different strain rates, 
which implies that some results provided in |3ZI are extrapolations of experimental data. 
The steady elongational viscosity rje is plotted versus strain-rate intensity Di in Figure 9 
(unfilled circles) together with its prediction by Eq. (82) (filled circles). For comparison, 
we also present the estimate of the elongational viscosity for a Newtonian fluid r]e = 3rjs 
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(asterisks). Figure 9 shows that Eq. (82) provides a reasonable correction of the Newtonian 
formula for the dependence of the elongational viscosity on strain rate. 

To further examine Eq. (82), we consider experimental data for dilute solutions of poly- 
acrylamide in water /glycerol solvent The observations depicted in Figures 4 and 6 of 
jin] show that the steady shear viscosity r)s and the first normal stress coefficient remain 
practically constant at shear rates up to 200-500 s^^ (the low-shear plateau). Formula (82) 
predicts that in this region of strain rates, the steady elongational viscosity rje should increase 
with strain-rate intensity (strain-thickening). This conclusion is in agreement with the data 
depicted in Figure 10 of |43j, which show a pronounced increase (by an order of magnitude) 
in the elongational viscosity with e in the interval of strain rates from 1 to 200 s~^. We do 
not provide qualitative comparison of the model predictions with these observations (as it 
was mentioned before, the assumption r = is not necessary fulfilled for polymer solutions, 
whereas an analog of Eq. (82) for r > becomes too cumbersome). However, it is worth 
noting that constitutive equations (42) can describe (at least, qualitatively) strain-thickening 
in polymer fiuids. This conclusion is rather surprising, because the strain-thickening phe- 
nomenon is conventionally associated with coil-stretch transition in polymer solutions, while 
the latter effect is not accounted for by the model explicitly. 

3. Our aim now is to show that Eqs. (42) can correctly describe the experimental data 
on elongational fiows with constant strain rates (in particular, strain-hardening observed in 
uniaxial elongational tests) when their parameters are found by matching observations in 
shear tests. For this purpose, we study the evolution of the transient viscosity 1]^ with time 
for the same LDPE melt that was used in the analysis of stress overshoot in shear tests. The 
experimental data in start-up elongational tests at Gq = 150 °C reported in j32] sue depicted 
in Figure 10. 

It is worth noting a pronounced difference between the observations in extensional and 
shear tests presented in Figures 8 and 10. While the stress overshoot and subsequent strain 
softening in shear tests (Figure 8) occur at times of the order of 1 to 10 s (at the chosen shear 
rates), a pronounced strain hardening is observed in extensional tests at times of the order of 
10~^ to 10~^ s. In this time interval, our assumption regarding rapid tuning of the number 
of strands to its steady value corresponding to a given strain rate becomes inadequate, and 
the evolution of the function N{t) should be taken into account. 

Within the reptation concept, the quantity N{t) equals the concentration of tubes sur- 
vived until instant t > 0. It was demonstrated in jSH] that changes in the number of survived 
tubes in a matrix of fixed obstacles are correctly described by a single-exponential function. 
This implies that, as a first approximation, the first-order kinetic equation can be introduced 
for the function N{t), 

where A^°(Z?i) is the steady concentration of strands at deformation with the strain-rate 
intensity D^, and Tt.{D{) is the characteristic time for changes in the number of strands at 
transition from the rest to the steady fiow. As the elastic modulus G is proportional to A^, 
we find from this equation that the function G{t) obeys the relation 

flG 1 
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where G°(A) is given by Eq. (43). 

After this correction of the governing equations driven by the necessity to account for 
the evolution of the elastic modulus, the transient viscosity rj^ in a start-up elongational 
test with a constant strain rate e is determined by Eq. (58), where the stress difference AE 
obeys Eq. (54), the function K,{t) is governed by Eq. (53), and the coefficient G is described 
by Eq. (83). 

The following procedure is applied to fit the observations at various strain rates e depicted 
in Figure 10. First, we set r = and calculate G^{e) and A{e) by using Eqs. (43) and (79) 
[where Di is given by Eq. (52)] with the parameters Go, Kq, i^g, Ka and va found in the 
approximation of observations in shear tests. As the initial condition Gin for Eq. (83) is 
unknown [it corresponds to the elastic modulus after sudden application of the external load 
and does not necessary coincide with Go in Eq. (43)], we begin with matching observations 
for the elongational viscosity in a test with the smallest strain rate e = 0.3 s~^ (for this 
test, the largest number of experimental data are presented in Figure 10). The curve r]^{t) 
is determined by 3 constants: Gin, & and Tr (formally, the coefficient h found in the shear 
tests can be utilized, but we chose to treat h as an adjustable parameter in order to compare 
appropriate results in extensional and shear tests: the closeness of values of h found in two 
different tests will confirm that the fitting procedure is stable). The quantities Gin, ^ and 
are determined by the same algorithm that was employed in Section 7. We fix some intervals, 
where the "best-fit" parameters are located, divide these intervals into subintervals, for each 
triad chosen from the subintervals integrate Eqs. (53) and (83) numerically (by the Runge- 
Kutta method with the time-step At = 1.0 • 10~^ s), and calculate r]^ from Eqs. (54) and 
(58). The best-fit values of Gin, ^ and Tr are chosen from the condition of minimum of the 
cost function ^ 

R-Y\vT''{tm)-rir\tm)\ , 

where the sum is calculated over all experimental points tm depicted in Figure 15, rjf^^ is 
the elongational viscosity measured in the test, and ry"""^ is given by Eq. (58). 

The initial modulus Gin ~ 2.0-10^ Pa ensures the best approximation of the experimental 
data in the extensional test with the strain rate e = 0.3 s~^. We now fix this value and fit 
other curves depicted in Figure 10 by using the same algorithm with only two adjustable 
parameters, b and Tr. Each set of observations is matched separately. 

Figure 10 demonstrates quite reasonable agreement between the experimental data and 
the results of numerical simulation at all strain rates under consideration. The adjustable 
parameter b is plotted versus strain-rate intensity in Figure 11. Comparison of the 
dependencies &( A) found by matching observations in shear (unfilled circles) and extensional 
(filled circles) tests reveals that the experimental values of b practically coincide, which 
implies that the fitting algorithm is stable. 

The parameter t,- is plotted versus strain-rate intensity in Figure 12 (unfilled circles). 
According to this figure, the parameter r,- strongly decreases with strain-rate intensity Di. 
The experimental data are approximated by the phenomenological relation 

log Tr = Tro - T,i log D;, (84) 

where the coefficients t^q and Tri are calculated by the least-squares method. Figure 12 
shows that Eq. (84) ensures fair approximation of the observations. For comparison, we also 
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present in Figure 12 the instants Tp, when the stress S12 in shear tests reaches its maximal 
values (filled circles). The observations are fitted by the formula similar to Eq. (84), 

log Tp = Tpo - Tpi log A (85) 

with the coefficients t^q and Tpi found by the least-squares technique. Figure 12 demonstrates 
that at all strain rates under consideration, the peak instants Xp substantially (at least, by 
an order of magnitude) exceed the relaxation times Tp. This confirms our assumption that 
the evolution of the elastic modulus G with time can be disregarded in the analysis of stress 
overshoot under shear. 

4. As the objective of this work is to study the rate-dependent behavior of polymer fiuids 
at large deformations, we do not concentrate on their response in oscillatory tests with small 
strains and relaxation tests. It is instructive, however, to make some comments regarding 
the model predictions in conventional viscoelastic tests. 

We begin with the analysis of governing equations at small strains. For this purpose, we 

set 

Be = I + 2ee, (86) 

where e is the Finger strain tensor for elastic deformation, substitute expression (86) into 
Eqs. (42), and neglect terms of the second order of smallness with respect to 4 and e, 
where e is the strain tensor for macro-deformation. Using Eqs. (43) to (47), we arrive at the 
expression for the stress tensor, 

5](t) = -P(t)I + 2Goee(t), (87) 

where the same notation P is employed for the unknown pressure. The evolution of the 
elastic strain tensor is described by the differential equation 

§W + -ee(t) = |(t), n = ^. (88) 
at Ts at Aq 

Formulas (87) and (88) reveal that Eqs. (42) are transformed into the Maxwell model 
at small strains, which implies that these relations can adequately describe observations 
in conventional oscillatory tests, provided that the one-mode model (42) is replaced by its 
multi-mode analog. As this procedure is straightforward, we do not dwell on it in the present 
study. 

Two comments are noteworthy. First, the fact that the elastic modulus at small strains 
coincides with Gq provides an additional opportunity to verify the model predictions, because 
Go can be measured independently in shear relaxation tests, on the one hand, and it can be 
calculated from Eq. (43) by fitting the experimental data on G{k) in start-up shear tests, 
on the other. Although we do not expect this procedure to result in the same values of Go 
for LDPE melt [we have only 4 data points for G{k) that are approximated by the three- 
parameter formula (43), which means that the accuracy of determining Go is rather poor], 
we would mention that the value Gq^"^ ""'^ = 6.5 ■ 10"^ Pa reported in Table 2 overestimates 
by less than twice the value Gq*^^ = 3.4 ■ 10^ Pa evaluated from the data in relaxation tests 
with small strains pp. 

Secondly, the relaxation time Tg in Eq. (88) is inversely proportional to the parameter 
^0- To derive constitutive equations that are equivalent to the Maxwell model at small 
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strains was the main reason to preserve the first term in phenomenological equation (46) 
despite the fact that all observations in the interval of strain rates under consideration can 
be approximated quite well by Eq. (79) without Aq. Following Eq. (88), we can introduce 
an effective relaxation time Ts for an arbitrary strain-rate intensity Gi as the reciprocate of 
the function A{Dj). The experimental data for Ts{Dj) are depicted in Figure 12 together 
with their approximation by the equation analogous to Eqs. (84) and (85), 

log Ts = Tso - Tsi log Di. (89) 

The coefficients Tsq and in Eq. (89) are found by the least-squares technique. Figure 12 
shows that the slopes of the curves Tj.{Di) and Ts(Di) are quite similar, and for all strain 
rates, the effective relaxation time Tg exceeds the relaxation time for the elastic modulus 
Tr by two orders of magnitude. This provides another confirmation of our hypothesis that 
the evolution of the elastic modulus G to its ultimate value G^{Di) may be neglected in the 
study of start-up tests. Figure 12 also shows that given a strain rate, the effective relaxation 
time Ts is substantially higher than the time for stress overshoot Tp, which means that Tp 
cannot be employed to assess the function A{D]), and the entire curves Si2(t) should be 
approximated to determined the model parameters. 

Finally, we consider relaxation tests with finite deformations. In these experiments, 
constitutive equations (42) with necessary correction (83) demonstrate a time-dependent 
behavior with two characteristic time-scales: (i) rapid relaxation driven by changes in the 
modulus G and (ii) slow decrease in stress associated with the Maxwell-type evolution of 
the elastic Cauchy-Green tensor Be in Eqs. (42). This implies that the stress in a step- 
strain test cannot be factorable into a product of a strain-dependent and time-dependent 
functions. As the same features are also revealed by the pom-pom model, we do not discuss 
them in detail, referring to |3T], where recent observations data are analyzed that confirm 
the presence of two time-scales in the viscoelastic response of polymer fluids. 

10 Concluding remarks 

A constitutive model has been developed for a polymer fluid. A polymer is thought of 
as an incompressible network of strands bridged by junctions. The network is treated as 
permanent (strands cannot separate from their nodes). The junctions between strands can 
slide with respect to their reference positions under deformation. No restrictions are imposed 
on the rate-of-strain tensor for sliding of junctions, whereas the vorticity tensor for sliding 
is proportional to that for macro-deformation. Stress-strain relations are derived by using 
the laws of thermodynamics. For a Gaussian network of strands, the constitutive equations 
involve three adjustable functions, for which phenomenological relations are introduced. 

In a broad sense, the constitutive equations may be treated as a combination of the 
Johnson-Segalman and Leonov models. The important advantages of these relations are 
that (i) they are based on the conventional multiplicative decomposition of the deformation 
gradient (no slip tensors are used), and (ii) they presume (in agreement with the classical 
theories in polymer physics) the strain energy of strands to be Gaussian (no complicated 
formulas are introduced for the elastic energy). 
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The governing equations are simplified for uniaxial extension and simple shear with finite 
strains. Explicit formulas are developed for the steady elongational and shear viscosities, as 
well as for the normal stress functions. 

To verify the model, three sets of experimental data are approximated for monodisperse 
polystyrene solutions with various molecular weights and molecular weights between entan- 
glements. Good agreement is demonstrated between the observations for stress overshoot 
in start-up shear tests and the results of numerical simulation at various strain rates. It is 
revealed that the phenomenological relations correctly characterize the effect of strain rate 
on the material parameters, and their coefficients change consistently with mass-average 
molecular weight and concentration of entanglements. 

To validate the governing equations, the evolution of stresses is studied in start-up exten- 
sional and shear flows on low-density polyethylene melt. It is demonstrated that constitutive 
equations (42) and (83) adequately describe changes in the elongation viscosity with time 
when their parameters are found by matching observations in shear tests. 

Our approximation of observations in extensional and shear flows of LDPE melt demon- 
strates higher quality of fitting compared to the pom-pom model. An advantage of consti- 
tutive equations (42) is that they can also be used to match the experimental data on dilute 
polymer solutions, a class of polymer fluids to which the pom-pom model is inapplicable. 

New formula (82) has been derived that expresses the steady elongational viscosity in 
terms of the steady shear viscosity and the first normal stress function. The vahdity of this 
equation is confirmed by observations on polypropylene melt. 

The following characteristic features of the constitutive equations have been revealed: 

1. The model can quantitatively describe overshoots for the shear stress and the first 
normal stress difference in polymer melts and solutions, as well as stress undershoot 
in polymer solutions. 

2. It can quantitatively describe strain hardening in start-up extensional tests. 

3. The governing equations can qualitatively predict strain-thickening in steady exten- 
sional tests. 

4. The model can qualitatively describe two-scale stress relaxation in polymer fluids in 
step-strain tests. 

5. At small strains, the constitutive equations are transformed into the conventional 
Maxwell model. 
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Table 1: Concentrations c, mass-average molecular weights M„, polydispersity indices -D, 
and average numbers of entanglements per chain M^/M^ for PS solutions and LDPE melt 



Abbreviation 




c 


M„g 


;/mol 


D 




Reference 


PSl 


3.0 


wt-% 


8.42 


■ 10*^ 


1.17 


10.0 


nn 


PS2 


0.1 


g/cw? 


1.09 


■ 10*^ 




6.1 




PS3 


0.1 


g/cw? 


8.24 


■ lO'' 




1.8 


M 


LDPE 






2.47 


■ 10^ 


15.2 


over 60 


MM 



Table 2: Adjustable parameters for LDPE melt and PS solutions 



Abbreviation 


logKA 




logK, 




Go Pa 


Kg 




r 


LDPE 


-1.09 


0.40 


-0.06 


0.0 


6.5 ■ 10^ 


0.172 


0.11 


0.0 


PSl 


-0.77 


0.40 


-0.75 


0.0 


79.7 


0.058 


0.34 


0.1 


PS2 


-0.37 


0.31 


-0.78 


0.60 


59.6 


0.144 


0.51 


0.2 


PS3 


-0.15 


0.37 


-0.32 


0.15 


22.0 


0.082 


0.91 


0.4 
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Figure 8: The shear stress S12 versus time t. Symbols: experimental data on LDPE melt 
in shear tests with the strain rates k = 10.0 (unfilled circles), 5.0 (filled circles), 2.0 
(asterisks), and 1.0 s~^ (stars) IJ. Solid lines: results of numerical simulation. 

Figure 9: The steady elongational viscosity rje versus strain-rate intensity Di. Unfilled 
circles: experimental data on polypropylene melt in uniaxial extensional tests [TFj . 
Filled circles: predictions of Eq. (82). Asterisks: the Newtonian estimate rje = 3rjs. 

Figure 10: The transient elongational viscosity rj^ versus time t. Symbols: experimental 
data on LDPE melt in uniaxial extensional tests with the strain rates e = 0.3 (unfilled 
circles), 1.0 (filled circles), 3.0 (asterisks), 10.0 (stars) and 30 s"^ (diamonds) |12]. 
Solid lines: results of numerical simulation. 

Figure 11: The dimensionless parameter b versus strain-rate intensity D^. Symbols: treat- 
ment of observations on LDPE melt in shear (unfilled circles) and extensional (filled 
circles) tests. Solid line: approximation of the experimental data by the constant 
6 = 0.85. 
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Figure 12: The relaxation time for the elastic modulus t — Tj. (unfilled circles), the time 
for stress overshoot r = Tp (filled circles) and the effective relaxation time for sliding 
of junctions r = Tg (asterisks) versus strain-rate intensity Di. Symbols: treatment of 
observations on LDPE melt in shear and extensional tests. Solid lines: approximation 
of the experimental data by Eqs. (84), (85) and (89) with t^q — —0.88, r^i = 0.48, 
Tpo = 0.47, Tpi = 0.62, and t^q = 0.40, r^i = 1.09. 
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